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We present a theoretical study of spin polarized transport in Fe atomic contacts using a self- 
consistent tight-binding Hamiltonian in a non-orthogonal s, p and d basis set, the spin-polarization 
being obtained from a non-collinear Stoner-like model and the transmission probability from the 
Fisher-Lee formula. The behaviour of an infinite perfect Fe wire is compared with that of an infinite 
chain presenting geometric defects or magnetic walls and with that of a finite chain connected 
to infinite one- dimensional or three-dimensional leads. In the presence of defects or contacts the 
transmission probability of d electrons is much more affected than that of s electrons, in particular, 
contact effects may suppress some transmission channels. It is shown that the behaviour of an infinite 
wire is never obtained even in the limit of long chains connected to electrodes. The introduction 
of the spin-orbit coupling term in the Hamiltonian enables us to calculate the anisotropy of the 
magneto-resistance. Finally whereas the variation of the magneto-resistance as a function of the 
magnetization direction is step-like for an infinite wire, it becomes smooth in the presence of defects 
or contacts. 

PACS numbers: 72.25.Ba,73.63.Rt, 75.47.Jn 



I. INTRODUCTION 

The electronic transport properties of atomic point 
contacts between two metallic electrodes have recently 
been the subject of great interest, both from the exper- 
imental and theoretical points of view. On the experi- 
mental side such atomic sized conductors can be obtained 
either by means of scanning tunneling microscopes [l[ or 
by using mechanically controllable break junctions Q. 
In some materials, for instance gold, monatomic chains, 
several atom long, are spontaneously obtained [l], Q . 

In these experiments the conductance is measured dur- 
ing the elongation process and shows clear plateaus cor- 
responding to the various stable atomic arrangements in 
the constriction region. In noble and alkali metals the 
last plateaus before breaking are very close to multiples 
of 2e 2 /h @| . Oscillations of the conductance with the 
length of the suspended chain are also observed [J- In 
transition metals even though plateaus are still present, 
the quantized behaviour is much less clear. However, 
thanks to the great mechanical stability of the break 
junction technique, the system can be stabilized for a 
given value of the conductance. This is essential in the 
case of magnetic materials since it allows to perform a 
series of magneto-resistive measurements revealing a re- 
markably large anisotropic magnetoresistance (AMR) ef- 
fects in materials like iron, cobalt or nickel 0, 0] . 

Understanding these transport properties is a challeng- 
ing problem and is of prime importance for their future 
applications in nano and spin electronics. Theoretically 
this relies on the calculation of the transmission factor 
T(E) of an electron at energy E from an electrode to the 
other in the ballistic regime which is relevant for such 
devices. In a perfectly periodic nanowire the transmis- 
sion factor of each electronic state is unity. Therefore at 



small bias voltage the conductance is determined by the 
number of linearly independent electronic states at the 
Fermi level propagating in a given direction (transmis- 
sion channels), and is quantized in units of 2e 2 /h or e 2 /h 
for a non-magnetic or magnetic material, respectively [7[. 

In a realistic system the electrodes are macroscopic and 
connected by an atomic sized constriction region. For 
such a system T(E) is limited by the number of channels 
in the narrowest part of the constriction. However, the 
transmission factor is usually smaller than this number 
due to scattering effects at the edges of the system. Sev- 
eral methods have been proposed to calculate the elec- 
tronic transport in atomic contacts but the most popular 
one is using an expansion of the electronic states on a lo- 
cal basis set which allows an easy partition of the system 
into three parts: the two leads and a central part. In this 
approach T(E) can conveniently be obtained from the 
Green function of the central part interactingwith the 
electrodes, by using the Fisher-Lee formula [8l Ifl Hoi. [lJ] . 

In this work we use a tight-binding Green function 
formalism to investigate the role of (geometric or mag- 
netic) defects and contacts on the transmission factor, 
and thus the conductance at low-bias voltages, of atomic 
junctions with simple geometries. Most of the calcula- 
tions we present here have been devoted to iron, using 
a realistic tight-binding Hamiltonian [l2| which includes 
magnetism and, possibly, the spin-orbit coupling term 
[l| ■ For the sake of comparison we have also considered 
gold which contrary to iron is non-magnetic and has elec- 
tronic states of pure s character at the Fermi level. 

We show in the following that the monatomic perfect 
wire model [|| [lj|, |T^, [l6| which has been used to explain 
experimental results is not reliable, at least for transi- 
tion metals, since the presence of defects and contacts 
strongly perturbs the transmission factor of the various 
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conductance channels. Indeed, even in the limit of long 
chains, the behaviour of the perfect infinite wire is not 
recovered due to contact effects. 

In Sec. II we present the theoretical model and give 
some details about practical computations. Sec. Ill A 
is devoted to the study of the influence of a geometrical 
or magnetic defect in a monatomic iron wire. The role 
played by atomic contacts is discussed in Sec. Ill B, and 
a comparison is made between Fe and Au. Conclusions 
are drawn in Sec. IV. Finally details of algebraic cal- 
culations are presented in the appendices. In appendix 
A we derive the Fisher-Lee formula for the general case 
in which the atomic spin-orbitals (and thus the overlap 
integrals) and the Hamiltonian matrix may be complex. 
Appendix B explains the method used for the calcula- 
tion of the surface Green function. Appendix C gives the 
main steps in the derivation of the transmission factor 
for a simple model studied in the text. 

II. METHOD 
A. The tight-binding model 

The electronic structure is derived from an Hamilto- 
nian operator H expressed in a non-orthogonal basis set 
of s, p and d atomic spin-orbitals hereafter denoted as 
\iXa) where i is an atomic site, A an atomic orbital and 
a the spin. The matrix elements (iXa\H\j/j,a') of the op- 
erator H in the \iXcr) basis form the matrix H. This 
matrix is written as the sum of four terms: 

H = Htb + -f^cc + -Hlcn + H so , (1) 

-Htb is a tight-binding Hamiltonian expressed in terms 
of Slater-Koster parameters which arc parametrized 
by a fit on ab-initio calculations for the bulk in 
the non-magnetic state [l2|. H cc accounts for the 
spin-polarization and is written as a simple intra- 
atomic Stoner exchange potential applied on d electrons 
— 1/2 Mtd-cr, where Mid is the d net spin magnetic 
moment vector of site i and u the Pauli matrix vec- 
tor. Note that this exchange potential allows for a 
non collinearity of spins. / is the Stoner parameter 
which we take equal to leV as in our previous studies 
on iron [HI, [T^j- -Hlcn is added to ensure quasi lo- 
cal charge neutrality in systems with inequivalent atoms 
[l3| . Its matrix elements (iAcr|i?LCN|j M 17 ') are written 
AP° n (A 9i + Aq^Sifdaa*, A pon is the penalization factor 
(in practice A pcn = 2.5eV), Aqi is the deviation of the 
Mulliken charge of site i from the valence charge, and 
S^Sau' arc the matrix elements of the overlap matrix 
S. Finally H so is the spin-orbit coupling (SOC) Hamil- 
tonian written as £,L.S where L and <S = fur/2 are the 
orbital and spin momentum operators, respectively. Only 
intra-atomic matrix elements between d spin orbitals are 
retained and £ is the spin-orbit coupling constant. It has 



been determined for Fc in rcf [l3| and is equal to 0.06cV 
when L and S are expressed in Bohr magnetons. 

The electronic (and magnetic) structure of a given sys- 
tem and its corresponding Hamiltonian are obtained by 
a self-consistent procedure since H cc and .Hlcn depend 
on the charges. 

This model has been checked on bulk, surfaces and 
monatomic wire of magnetic iron, and good agreement 
with ab initio calculations [l3| was found. Finally let us 
mention that, although we have shown in recent publica- 
tions that orbital polarization effects EH (taken into 
account by a more accurate expression of the electronic 
interaction Hamiltonian H cc ) may have an influence on 
the electronic and magnetic properties of nanostructures, 
we have ignored them in the following since the aim of 
the present work is to discuss the effect of contacts on 
transport properties that are present whatever H cc . 

In the following, for the purpose of deriving qualitative 
arguments, we have also used an Hamiltonian H limited 
to Htb for an s band. 



B. Electronic transport formalism 

We calculate the electric conductance in the Landauer- 
Biittiker approach 0- 
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FIG. 1: Schematic atomic model for the transmission calcu- 
lation in the Landauer approach. A scattering region (C) is 
coupled to two semi-infinite leads (L and R), L and R are 
assumed to interact only through C. 

The system is divided into three parts (Fig. [I]): a 
central scattering region (C) and two semi-infinite leads 
(L and R) with two-dimensional periodicity. 

The transmission probability T(E) of an electron with 
energy E is obtained from the Green function formalism. 
Our approach is close to the one proposed by Sanvito et 
al. [Kj. The scattering region in Fig. [T]is chosen wide 
enough so that no direct interactions remain between the 
two leads. The complete Hamiltonian and overlap ma- 
trices of our system are then 
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(2) 

The matrix He is of size iVc x Nq where Nq is the 
number of atomic spin-orbitals in the central region. ifj, 
and Hr are semi-infinite. The number of non-zero ele- 
ments of the coupling and overlap matrices H a c and S a c 
is N a x Nc where N a is the number of spin-orbitals on 
atoms in lead a (L or R) that have hopping and overlap 
integrals with the central region. 

Here we have considered the general case where the 
overlap matrix (thus the atomic spin-orbitals), and the 
Hamiltonian matrix elements may be complex. The ex- 
pression of T(E) is derived in appendix A following the 
same spirit as in the work of Viljas et ah In this ap- 
pendix we show that the Fisher-Lee relation is still valid 
in the general case, i.e.: 



T(E) = TiiT L (E)G c (E)T R (E)Gl(E)), 



(3) 



where Gc{E) is the retarded Green function of the cen- 
tral part defined by the relation: 
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where I a is the identity matrix and E + = lining E + irj. 
In practice we add a small but finite imaginary part r\ 
to the energy. In the following the superscript + will be 
omitted for convenience. 

Solving Eq. I for G C (E) yields: 

G C (E) = (ESc-H c -'Z L (E)--ZR(E)y 1 (5) 
T a (E) = i(E a (£0 - St (25)) (6) 

Si and are the self energy terms which account for 
the coupling of the central part to the leads, i.e.: 



The calculation of the central Green function Gc is now 
reduced to the inversion of a matrix of size Nc x Nc- 

Finally, experiments being performed at small bias 
voltage and low temperature [y] Eq |A22l leads to a con- 
ductance given by 



g = -T(E F ) 
h 



(8) 



where Ep is the Fermi level of the system. Note that the 
quantum of conductance is here given by e 2 /h and not 
2e 2 /h because we are considering magnetic materials. 



C. Computational details 

The atomic structure of the leads is that of a semi- 
infinite crystal with two dimensional periodicity but, due 
to the presence of the central region, their electronic 
structure and magnetic moments are modified near the 
scatterer and no longer periodic. However the choice of 
the scattering region is arbitrary and can include a part 
of the leads containing the atomic planes around which 
the matrix elements of the Hamiltonian are significantly 
different from the bulk ones. We therefore define an ef- 
fective scattering region by adding a part of the leads to 
the central part. In practice for bec (001) three atomic 
planes are added on the left-hand side of the scatterer 
and two on the right-hand side. The computing process 
runs as follows. First we calculate the electronic structure 
and magnetic moments of the effective scattering region, 
and since in most of the cases considered in this work the 
left and right leads are identical, we have added periodic 
boundary conditions in the three directions. With that 
procedure the bulk like behavior in the additional atomic 
planes is recovered faster. We thus obtain the Hamilto- 
nian of the central region He from the self-consistent 
calculation in the effective scattering region, the Fermi 
level being fixed to that of the leads. 



Bulk Bulk Bulk S-N 



Left lead 



S-N Bulk Bulk Bulk 



Effective scullering region 



Right lead 



FIG. 2: Definition of the scattering region for practical cal- 
culations. 



H a (E) = (ESl c - Vl c )g s a (E)(ES aC - V aC ) (7) 

where g^{E) is the surface Green function of the uncou- 
pled lead a. This surface Green function is calculated in 
the absence of the scattering region and thus is different 
from the Green function G a defined in Eq. 2J The func- 
tions ga(E) need to be calculated only on the atoms of 
the leads which are in contact with the scattering region. 



In this geometrical configuration g^{E) (in Eql7|) is 
the surface Green function of a truncated bulk in which 
the layer Hq and intcrlaycr Hi matrix elements (See Ap- 
pendix B Eq. IB1I ) of the Hamiltonian are the same as 
in the bulk, i.e., they are not modified by surface self- 
consistency effects. The Hamiltonians of the leads H^ 
and Hr are obtained from a bulk calculation. We then 
calculate the transmission T(E) (Eq. [3]). 
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III. RESULTS 
A. A test case: the monatomic iron wire 

We have first considered transport in the simple test 
case of a monatomic wire along the z axis where the re- 
sults can easily be interpreted. Moreover, this case is rele- 
vant for the study of break junctions. Indeed experiments 
on gold or platinum have shown that a short monatomic 
wire may appear and stabilize during the breaking pro- 
cess d, [l9[. In iron junctions, both theory and exper- 
iments do note show evidence for the appearance of a 
wire [2(j ■ Nevertheless the contact is still atomic and the 
geometric environment of the contact atom is very close 
to the environment of an atom in a monatomic wire. 

In the following we study the influence of various de- 
fects on the transmission probability T(E). 



1. Geometric defect 

First we study the effect of a geometric defect on the 
transport properties by considering a perfect infinite wire 
in which the distance between two atoms is stretched and 
SOC is neglected. The leads arc two semi-infinite wires 
with an interatomic distance do = 2.27 A and magnetic 
moment mo = 3.21^b at equilibrium. The effective scat- 
tering region is made of 10 atoms at the same equilibrium 
distance, save for the distance d between the 5th and the 
6th atoms which is stretched from equilibrium to 3.25A 
(Fig. [3]). The magnetic moments in the scattering region 
are computed sclf-consistcntly. Our calculation showed 
that the magnetic moments of the atoms near the defect 
increase with the distance d as expected. 




E-E p (eV) 



FIG. 3: Transmission T(E) for an iron monatomic wire with 
the geometric defect shown in the upper part (d = 2.27A). 



The electronic transmission for different stretchings is 
reported in Fig. [3J For the perfect wire (black curve), 
each channel has a transmission probability equal to 
unity and the total transmission is simply given by the 
number of electronic states at energy E. As the wire is 
stretched, the transmission of each channel decreases. 

We can see that the reduction of transmission proba- 
bility is strongly energy dependent. The channels corre- 
sponding to d xy and d x 2_ y 2 orbitals (5 bands), which are 
localized in the plane perpendicular to the transmission 
direction z, are more affected by the stretching of the cen- 
tral atom than the channels corresponding to s and d z 2 
orbitals. This explains the flattening of the high trans- 
mission peaks at —3eV and 0.25eV which correspond to 
S orbitals. 



2. Magnetic defect 



We now look at the effect of a domain wall on the 
electron transport in the absence of SOC. 




FIG. 4: Transmission T(E) of an iron monatomic wire in the 
presence of a N atom long magnetic wall (see text). 

The leads are two semi-infinite wires with opposite 
magnetization directions. The scattering region contains 
a magnetic wall of length N. The angle between the 
magnetization directions of two consecutive atoms in the 
wall is chosen to be constant and equal to tt/ (N + 1). 

A first calculation shows that the presence of an abrupt 
wall (N = 0) results in an important loss of transmission. 
The transmission T(E) = 2 is due to the s orbitals which 
are the only channels able to transmit in this case. Indeed 
in transition metals, the magnetic moment is mostly car- 
ried by the d electrons and in an iron wire the d orbitals 
are totally polarized. Thus no transmission is possible 
through the wall. On the other hand, s orbitals are only 
partially polarized and are able to transmit. 

When the magnetization is flipped by 7r/2 on one atom 
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between the leads, leading to a domain wall with N = 1, 
the d channels open and the transmission increases. This 
spin-flip transmission is made possible by the states in the 
flipped atom that mix spin up and spin down orbitals. As 
the length of the wall increases, the transmission quickly 
tends towards the transmission of a perfect wire with no 
domain wall. 

Thus, in the ballistic regime, a magnetic wall a few 
atom wide has almost no effect on the electron trans- 
mission. This results is in good agreement with the cal- 
culation of Vclev and Butler [2l[ who showed that large 
magneto-resistance is expected only in contacts with very 
narrow domain walls. 



transmission is no more equal to unity. The curve of the 
transmission as a function of the magnetization direction 
for the same type of defect as in Sect MI ATI is no longer a 
step between two integer values but becomes continuous 
with a maximum reached at 9 = tt/2 (Fig. [5]). The 
AMR is softened by the presence of the defect as already 
reported by Jacob et al. [22| ■ We have seen in Sec. IIII A II 
that the 5 channel was the most affected by a geometric 
defect in a wire. Since the AMR in an iron wire is only 
due to electrons in these S bands it is strongly affected 
by the presence of a defect (Note that the domain of 
stretching is different in Fig. [5] and in Fig. [3]) . 



3. Influence of spin orbit- coupling 
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FIG. 5: Transmission at the Fermi energy T(Ef) in an iron 
monatomic wire with the geometric defect shown in Fig. [3] 
as a function of the magnetization direction 8. 

In this section SOC is taken into account. Therefore 
the electronic structure depends on the magnetization 
direction and anisotropic effects appear. In a perfect 
monatomic wire of a transition metal, SOC removes the 
degeneracy of the S bands except when the magnetiza- 
tion is perpendicular to the wire (9 = tt/2) [l3j . In iron, 
these bands lie around the Fermi level and the splitting 
can be strong enough to bring one of the bands above it 
when 9 varies. Thus, the conductivity of the wire (which 
is proportional to the number of bands crossing its Fermi 
level) decreases by 1 as the magnetization orientation 
is switched from 7r/2 to (Fig. [3 black curve). This 
anisotropic magneto-resistance (AMR) in 3d transition 
metal wires has already been reported [H, [l5[ . Note that 
here and in all the following calculations in which 9 is 
varied, a penalization function has been added to the 
Hamiltonian in order to keep a fixed value of 9 during 
the self-consistent iteration process (See Ref. fl3(| ) . 

If a geometric defect is present in the wire, the channel 



B. Influence of the atomic contacts 

We have seen in Sec. IIII Al that the transmission of 
an iron atomic wire is affected by geometric or magnetic 
defects. Another strong effect arises since in practice the 
wire is always connected to some leads. These contact 
effects can become very important for instance in break 
junction experiments where the leads are bulk-like while 
the scattering region is quasi one-dimensional. To have 
a better understanding of the role of the contacts in the 
Landaucr formalism applied to break junctions, we stud- 
ied the transmission of some simple systems. 



1. The contact effect in a s band model 



L C R 
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b) 



FIG. 6: Geometries used in our calculations, a) N atom 
long wire with hopping integral j3 connected to two semi- 
infinite wires with hopping integral j3o by means of the contact 
hopping integral 7. b) N atom long wire between two fee (001) 
surfaces. The atoms 1 and N are in four-fold position relative 
to the surface. The distance between atoms in the wire, in 
the leads and at the contact is given by the bulk interatomic 
distance at equilibrium of the considered material. 
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The simplest way to model the contact is to consider 
a non-magnetic wire of N atoms with a single s orbital 
per site and an hopping integral f3 between first nearest 
neighbours. This wire is connected to two semi-infinite 
wires which have a different hopping integral /3q but the 
same atomic level. The contact is established by a hop- 
ping integral 7 between the leads and the edge of the 
finite wire (Fig. [5] a). The three hopping integrals are 
taken negative, overlap integrals are neglected, and all 
first neighbour interatomic distances are equal to unity. 
The transmission of this system can be calculated ana- 
lytically (see Appendix [C]) and is given by: 




FIG. 7: Transmission of the system shown in Fig. [6] a) 
calculated from Eq. [9] as a function of energy and of the 
number N of atoms in the central wire, with /3o = —2, /? = —1 
and 7 = —1.2. b) as a function of energy and of the strength 
I7I of the contact, with j3o = —2, f3 = — 1 and N = 5. 

T N (E) = ( 9 ) 
2 sin k sin kn 

2 -5-J5 " 

^ sin(iV - l)k - 2 sin Nk + ^ sin(iV + l)fc 

where Tjv (E) is the transmission coefficient of a N atom 
long wire at energy E with E = 2(3 cos k = 2/3oCOsfco- 

In the following we study the case \{3\ < \(3o\, i.e., the 
dispersion of the central part is narrower than that of 
the leads. In Fig. [7^ we plot the transmission Tn(E) 
when \/3\ < \j\ < \(3q\, for different values of N. Note 
that k is imaginary when E > 2\(3\. In the limit of large 
N the transmission is non-vanishing only inside the en- 
ergy band of the N atom wire ( -2|/3| < E < 2|/3| ). 
Inside this band, the transmission oscillates and presents 
a number of peaks that increases with N (there are N 
peaks for a wire of N atoms for this specific choice of (3, 
(3q and 7). It is interesting to note that, in the middle 
of the band (where E = and k = ko = tt/2), the trans- 
mission oscillates between two values when N goes from 
odd to even. If N is odd, TWO) = 2 and if N is even 
Tjv(O) = 8( 7 2 /(/W + (PPo)/ 1 2 )' 2 - When 7 = ~VM, 
this oscillation disappears. Finally for N small (N < 5) 
the transmission has a significant exponentially decaying 
tail outside the energy range [— 2|/3|, 2\/3\]. 

In Fig. [7b, we plot the 7 dependence of the trans- 
mission of a five atom wire. For I7I = y/(3f3o, the trans- 
mission is very close to 2 on the whole energy spectrum. 



For small values of N sharp peaks appear at energies 
close to the eigenvalues of the isolated finite chain. When 
I7I increases the peaks broaden and are shifted so that 
some of them disappear. 

Thus, the number of oscillations in the transmission 
curve increase with the number of atoms in the wire and 
the strength of the contact 7 controls the amplitude of 
these oscillations. These results show that a too strong 
or too weak contact kills the transmission while with an 
appropriate choice of 7 the transmission can be perfect. 

Let us now consider a more complex geometry where 
the leads are two semi-infinite surfaces as shown in Fig. 
[BJd. We illustrate this case for gold using first a basis set 
limited to overlapping s orbitals with a parametrization 
taken from Mehl and Papaconstantopoulos [l2||. Note 
that however we assume that the onsite term does not 
vary with the atomic environment and is chosen as the 
zero of energies. This model can be used as a first ap- 
proach since the electronic conduction in gold is domi- 
nated by the s electrons as Ep belongs to the s band. 
The choice of gold is especially interesting because ex- 
periments on gold atomic contacts have shown its abil- 
ity to form monatomic wires several atom long. Fur- 
thermore the conductivity of these wires exhibits specific 
features such as quantization and parity oscillations [H[ . 
The inter-atomic distance in the wire and in the con- 
tacts is the nearest neighbour distance in bulk gold at 
equilibrium (2.88A), and the surfaces of the two leads 
are fcc(OOl). The results are shown in Fig. [H^i. As for 
the perfect infinite wire we observe that the transmission 
curves oscillate as a function of N and E. In the middle 
of the band, the transmission is close to 2 and exhibits 
even-odd oscillations. 

Additional calculations with different distances be- 
tween the finite wire and the surfaces revealed that, as in 
the infinite wire case, the strength of the contact affects 
the amplitude of the oscillations. From these results, we 
can conclude that in the case of a pure s band, the choice 
of wire leads or bulk-like leads does not much modify the 
transmission of the system. 



2. The contact effect in metals with spd bands 

To investigate the contact effect on a system with a re- 
alistic electronic structure we have calculated the trans- 
mission for a gold wire of N atoms connecting two (001) 
gold surfaces using the spd TB parametrization of Mehl 
and Papaconstantopoulos [l2j, and neglecting SOC. The 
results are shown in Fig. [SJd. We can see that in the 
energy range where only the s band is present, the trans- 
mission curves are similar to Fig. [8^, with oscillations 
as a function of N and E and a transmission close to 2. 
In the energy range corresponding to d states T(E) has 
a complex behaviour which does not show simple oscil- 
lations. It is clear that the presence of contacts strongly 
affects the transport of d electrons. However in gold the 
d electrons do not participate in the conductance at low 




o 
o 



o 



FIG. 8: Transmission T(E) of a N atom Au wire between two 
Au(OOl) surfaces, a) s parametrization. b) spd parametriza- 
tion. Note the complex behaviour in the d band energy range. 
Ef is the zero of energy. 



bias voltage since the Fermi level has a pure s charac- 
ter. Thus a gold atomic contact exhibits a quantization 
of its conductance and odd-even oscillations as observed 
experimentally 0] and explained theoretically [24| ■ 

On the opposite, in transition metals like iron, s, p 
and d electrons are present at the Fermi level and the 
behaviour of the atomic contact conductivity is more in- 
tricate 25]. To illustrate qualitatively the importance of 
contacts in such a case we have calculated the transmis- 
sion of a non-magnetic iron wire of 5 atoms connecting 
two bec (001) leads and decomposed it into the contribu- 
tions of the different transmission channels. The results 
are shown in Fig. [9J We have also reported the local den- 
sity of states (LDOS) on the central atom of the wire, on 
the surface of the leads when disconnected from the wire 
derived from g^{E) as explained in Sec. Ill CI and in the 
bulk. For the sake of comparison we have calculated the 
same quantities in the case of a perfect infinite wire. 

The <7, 7r and 5 decomposition corresponds to the three 
band symmetries found in the perfect wire and defines its 
transmission channels. The wire being along the z axis, 
the a channel results from the hybridization of s,p z and 
d z 2 orbitals, the n channel homp x , p y , d xz and d yz , and 
the 5 channel from d x z_ y 2 and d xy orbitals. 

In the case of the wire between two surfaces a, tt and 
5 orbitals become coupled. Accordingly the transmission 
channels have no longer a pure a, 7r or 5 character. The 
general way of determining the transmission channels is 
to calculate the t matrix defined as 
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FIG. 9: Bulk (a), local DOS at the surface of the leads (b) 
and on the central atom (c) of the wire, and transmission (d) 
in an infinite Fe wire (left) or in a 5 atom Fe wire connecting 
two bcc(001)Fe surfaces (right). The contributions of the a 
(black), 7r (red) and 5 (green) channels (see text) have been 
resolved. 



(10) 



such that T(E) = Tr(W). The transmission channels are 
the eigenvectors of tv and the correspon ding eigenvalues 
are the associated transmission factors [26|, |22|. Since 
the choice of the scattering region C is somewhat arbi- 
trary we took the central atom of the finite wire. The 
calculation shows that for a wire as short as 5 atoms the 
transmission channels have nearly a pure a, n or S char- 
acter. However the value of the transmission factor is 
no longer quantized but presents sharp oscillations as a 
function of E (See Fig. [9J. In order to obtain a physical 
insight into this result we also show in Fig. [9] the sum of 
the LDOS projected on the various orbitals contributing 
to each channel. For the central atom of the wire the 
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LDOS have some similarities with the ones of the perfect 
wire, especially the 6 LDOS. Indeed all the weight of the 
S LDOS is concentrated in the same energy range as in 
the perfect wire. This is no longer true for the LDOS of 
the leads, since the symmetries are completely changed. 
As a result the various LDOS have a significant weight on 
the total energy spectrum and are thus lowered. Conse- 
quently the transmission channels are strongly affected. 
In particular the contribution of the 6 channel is almost 
zero. 

Finally let us note that the connection between the 
lead and the wire considered here is very abrupt and a 
possible origin of the weak transmission could be this 
unrealistic geometry. However we do not believe in the 
pertinence of this suggestion since we have seen that for 
a pure s band this geometry can lead to almost perfect 
transmission. Moreover similar calculations for magnetic 
Fe in an spd model, where we assumed a smoother con- 
tact (small pyramids), gave also low channel transmission 

3. Anisotropic magnetoresistance effects in Fe atomic 
contacts 

We have seen in Scc lIII A3l that an infinite iron 
monatomic wire shows a step-like variation of T(Ep) as 
a function of the magnetization direction (Fig. [5]). This 
abrupt AMR takes its origin in the behaviour of the 8 
band and is softened by the stretching of one bond. In 
the previous Sec. ((Ill B 2p we have also demonstrated 
the influence of contacts on the transmission probability 
of each channel which may almost vanish. It is there- 
fore interesting to investigate the role of contacts on the 
AMR. 

We have calculated the transmission at the Fermi 
energy in an iron wire of N atoms connected to two 
bcc(OOl) Fe leads for different magnetization directions. 
The results are shown in FigfTHlas a function of N. 

The AMR effect is still present, but the curves are con- 
tinuous and the difference of transmission for magnetiza- 
tions along or perpendicular to the wire is lower than 1. 
The switching between two quantized conductance val- 
ues predicted in the perfect infinite wire (Fig. [5l black 
curve) disappears. 



IV. CONCLUSION 

The idealized model of an infinite monatomic wire has 
often been used to interpret the conductance measure- 
ments in atomic contacts. However, even though a wire 
may exist between two electrodes, it may present some 
defects. Furthermore this wire is necessarily finite and 
the contact between the two leads cannot be ignored. In 
this work we have investigated, on the one hand, the in- 
fluence of structural (stretching of a bond) and magnetic 
(domain wall) defects on the transport properties of an 




e 



FIG. 10: Transmission at the Fermi energy of a N atom long 
Fe wire connecting two bec (001)Fe surfaces as a function of 
the magnetization direction. 



infinite wire. On the other hand, we have studied the 
contact effect, both in a ID model and in a more elab- 
orated geometry in which the leads are represented by 
semi- infinite 3D crystals (See Fig. [5]). 

The role of defects in a Fe monatomic wire has been 
computed with an elaborate tight-binding Hamiltonian 
written in an spd basis set. As expected the stretch- 
ing of one bond or the presence of a magnetic domain 
wall decrease the transmission factor (which is no longer 
quantized) the more as the defect is pronounced (large 
stretching or abrupt wall) except outside the range of 
the d bands (i.e., between 2 and 4 eV above the Fermi 
level). Surprisingly the effect of a domain wall almost 
vanishes as soon as its width overcomes three atoms. 
When spin-orbit coupling is taken into account, the step- 
like behaviour of T(Ep) as a function of the magnetiza- 
tion direction 9 (i.e., the AMR effect) is softened by the 
stretching of one bond in a perfect infinite wire. 

The influence of contacts has been compared for two 
metals: Au and Fe in which the electronic states at the 
Fermi level have s and d character, respectively. In gold 
the effect of contacts on the transmission factor is strong 
in the energy range corresponding to d bands, but rela- 
tively limited around the Fermi level. Note however in 
the latter case the presence of smooth odd-even oscilla- 
tions. 

In iron the transmission factor presents very sharp os- 
cillations as a function of energy in the presence of con- 
tacts, and some transmission channels can even almost 
disappear, in particular the channel of d xy and d x i_ y i 
character (S band). This latter result casts some serious 
doubts on the interpretation of the AMR effects in mag- 
netic contacts based on the band structure of a perfect 
infinite wire. In conclusion we have demonstrated that 
contact effects strongly change the electronic transmis- 
sion of a monatomic wire in transition metals. Thus it 
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is hoped that by modifying the contacts one could tailor 
the transport properties of nanowires. 
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APPENDIX A: DERIVATION OF THE FISHER 
LEE FORMULA IN THE GENERAL CASE 

The aim of this appendix is to prove the Eql3] of the 
main text which gives the electronic transmission in the 
general case where the basis set is non-orthogonal and 
not necessarily real and in which the matrix elements of 
the hamiltonian may be complex. 

We consider an electron which at time t=— oo is de- 
scribed by an incident Bloch wave ijj^ of the left lead, 
i.e., obeying: 



(Hi - E {n) S L )i 







(Al) 



and switch on adiabatically the connection to the scat- 
tering region and the right lead. At a time t< the 
hamiltonian is thus H(t) = + (He + Hr) exp(ei), e is 
a small positive factor which ensures that the connection 
is established adiabatically and ^>(r, t) can be expanded 
as a linear combination of atomic spin-orbitals with time 
dependent coefficients a,( n \ a (t). At t=0 the wave func- 
tion ip(r,t) has evolved towards a stationary state of H 
with the same energy as the incident wave that will be 
denoted in the following as t/j for short. 

The charge Nq contained in a given part CI of the 
system (Cl = L(left lead), C(scattering region), i? (right 
lead), L + C, C + R) at time t is defined according to 
the Mulliken population analysis and is written as (after 
having summed over all possible incident waves): 



2 K dt dt ' 2 V dt dt 



since the overlap populations between L+C and R cancel 
in this difference. Let us now calculate -g^ 1 : 



(A5) 



+ Re 

«occ,an,/3f!/as2 

From the time dependent Schroedinger equation writ- 
ten in the non-orthogonal basis set we have: 



(A6) 



d 

nP Qj} a (n)p(t)}- 



Substituting Eq |A6l for 



c)a 



the right-hand side of Eq lA5l yields 



into the first term of 



dt 



d 

ihS anfin —}a (n) ^{t) + c.c. 



(A7) 



Cl + Cl denoting the whole system. At t — the state n is 
a stationary state of energy £"(„), thus daf n -\ a /dt can be 
replaced by — j-E( n )a( n ) a (0). Finally substituting Eq |A7l 
for dN'Jdt in Eq|A4]we find: 



N n (t) = Re[ £ (t)S{, O'oJ 

= ReTi n P{t)S (A2) 

where an is an atomic spin-orbital belonging to the part 
CI and (3 is any spin-orbital of the system, S is the over- 
lap matrix and P a p{t) = J2n occ a (n) a (*)«(„)/?(*) arc tnc 
elements of P. Note the presence of the real part in Eq. 
A2I which is necessary when the hamiltonian or the basis 
set is complex. In the following the P matrix is divided 
into blocks, similarly to H and S (see Eq|2]of the main 
text). It is convenient to define another population JVq 
which is obtained by retaining in Eq |A2l onlv the overlap 
populations between orbitals included in CI, i.e.: 



NUt)=TcP a (t)8 t 



(A3) 



The current Il^r, coming from the left is given either 
by the electrons flowing from L + C or by those entering 
in the right lead R. Thus: 



h^R = -^r c W CR (E {n) )P RC + c.c. (A8) 



with: 



W(E {n) )=H-E (n) S. 



(A9) 



The response of the system to iM„) is denoted as tp', 
tp = ip l i n ) V*' obeying the Schroedinger equation of the 
whole system: 



(H -E [n) SM n) +^') = 



(A10) 



In the non-orthogonal basis set the vectors ip l , n ~. and 
ip' can be decomposed into blocks corresponding to the 
three parts of the system (L, C, R), i.e., 




a (n)L 

and tp = | a( n ) C 

a (n)R 
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respectively. Substituting these components into Eq lAlOl 
with Hlr = ^lb. — an d taking Eq lAll into account 
yields after simple algebraic manipulations: 



a (n)c = Gc(E {n) )W C L(E {n] )^ n)L (All) 



a ( „)H = G R c(E {n) )WcL(E {n) )a\ n)L (A12) 



with: 



G(E {n) ) = (E [n) S-U)-\ (A13) 
Finally from EqfJJof the main text it is easy to show that: 



the square matrix of dimension Nl x Nl with ele- 



ments Y^„ a \n) aL a ln)p L S ( E ~ E (n)) which is nothing but 
i(gL(E) — g? L (E))/2n (with obvious notations) in the gen- 
eral case where H and the basis set may be complex. By 
making use of the cyclic invariance of the trace we obtain 
finally: 



1 r°° 

-J Tr(Gl(E)T R (E)G c (E)T L (E)) 



x f(E-E FL )dE 



(A20) 



where / is the Fermi function. 

The current due to an incident wave from the right lead 
Ir^l is obtained by interchanging the R and L indices. 
The total current is 



G RC (E {n) ) = gR (E in) )W RC (E (n) )Gc(E {n) ) (A14) 
with: 



1 = 1 



SR(E(n)) — (E( n )S R - Hr) . 



(A15) 



Substituting Eq |A14l for G RC into Eq |A12l and using 
Eqs lAllI and IA12I to obtain P R c leads to: 

= i^[a^ )i W £C (S (n) )G^(£; ( „ ) )W c7il (£; ( „ ) ) 

n occ 

x g R (E {n) )W R c(E {n) )G c (E {n) )W CL (£(„) )a{ n)£ 
- c.c] (A16) 

where G^ is the hermitian conjugate of G, or: 



ll^r 



^J2K)L W ^(E {n) )Gl(E {n) ) 



x W CR (E {n) )(g R (E {n) ) - g R (E {n) ))W RC (E {n) ) 
x G c (E {n) )W CL (E {n) )^ n)L }. (A17) 

From Eqs[7]and[6]of the main text we have: 

W CR (E {n) )( SR (E {n) ) - gl(E {n) ))W RC (E {n) ) = 

-iT R (E (n) ) (A18) 

thus: 



il->r 



■is 



^)l^lc(E)gUe) 



x r R (E)G c (E)W CL (E) 
x j {n)L 5{E - E {n) )AE 



(A19) 



where Ep L is the Fermi level of the left lead. The 
integrand can be written as a trace if we introduce 



L^R 



I 



R^L 



(A21) 



and, since the current should vanish whatever the system 
when Ep R = Ep L : 

I = ~ j T(E) x <J(E -E Fl )- f(E - E Fa ))dE. 

(A22) 

with 

T(E) = Tr(T L (E)G c (E)T R (E)Gl(E)) (A23) 

T(E) is the electronic transmission given by the Fisher- 
Lee relation (EqJ3]of the main text). 
Actually it can be shown that [28| : 



Tt{T l (E)Gc.{E)T r {E)G^{E)) = 
Tr(T R (E)G c (E)T L (E)Gl(E)). 



(A24) 



APPENDIX B: CALCULATION OF THE 
SURFACE GREEN FUNCTION 

To calculate the Green function at the surface of the 
lead ga(E) (See Sec III Bl and Sec. IHCj) . we use an iter- 
ative scheme. Each lead has two dimensional periodicity 
and can be viewed as a semi-infinite succession of iden- 
tical layers (Fig. ITT]) with the same interlayer spacing. 
The layers contain several atomic planes so that there 
are interaction only between consecutive layers. 

The Hamiltonian of the semi-infinite lead can thus be 
written as 



H a = 



\ 







H H\ 

H\ H q if x 
H\ Ho H x 
H ] Ho J 



(Bl) 
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FIG. 11: 

repeated. 



Each lead can be built from a layer periodically 



where Hq is the Hamiltonian of a layer and Hi is the hop- 
ping matrix between two successive layers. The overlap 
matrix S a has the same form. Hq and H± are matrices 
of size N a where N a is the number of atomic spin-orbitals 
in a layer of the lead a. 

We need to calculate the surface Green function only 
on the atoms that have interactions with the scattering 
region. Consequently we just have to calculate the sur- 
face Green function on the last layer. Thus is a matrix 
of size N a x N a defined by the relation: 



APPENDIX C: TRANSMISSION OF A FINITE 
MONOATOMIC WIRE CONNECTED TO TWO 
SEMI-INFINITE WIRES 



We consider two identical semi-infinite linear chains 
(leads) connected to a finite chain of N atoms of a dif- 
ferent species (scattering region) (See Fig. [BJi). All the 
atoms are non magnetic and equally spaced with an in- 
teratomic distance taken as unity. An s atomic orbital 
is centered on each atom. The electronic states are ob- 
tained from a tight-binding scheme in which the non- 
orthogonality of the s orbitals is neglected. Furthermore 
all the atomic levels of the system are assumed to be the 
same and taken as the energy zero. The hopping integrals 
are limited to first nearest neighbours and are equal to 
/3q and (3 in the leads and scattering region, respectively. 
Finally the connection between the leads and the scat- 
tering region is established by an hopping integral 7 (See 
Fig. [S^). The Schrocdingcr equation projected on each 
atomic site yields (for each spin): 



\ 



ESq — Hq ESi — Hi 

ES\ - H{ ESq - H Q ESi - Hi 
ES\ - H[ ES q -H J 



\ ( 



G 



s-i 



G 



s-is 



r<ss-i „s 



In„ 

i Na 

I Na J 
(B2) 



Now we add a layer (denoted as 5 + 1) to the sys- 
tem. Using the Dyson equation, we can calculate the 
surface Green function g^ +1 (E) of the new system from 
the surface Green function of the old system and the 
Hamiltonians Hq and Hi: 



,s+i 



(E) 



{ESq -H Q - (ES\ - H\)g s a (E)(ESi - Hi))' 



(B3) 



Due to the periodicity of the system, adding a layer 
does not change the Green function on the surface layer. 
Therefore, g^ +1 (E) = g^{E) and the surface Green func- 
tion is solution of the equation 



9 S a (E) 
(ES - 



Hq 



(ES\ - Hl)g s a (E)(ESi - Hi)) 



-1 (B4) 



This equation is solved iteratively using the quick it- 
erative scheme proposed by Lopez Sancho et al. [29[ . 



0oa n -i ~ Ea n + P a n+ i = 
A)fl-i — Eciq + 701 = 
7a - Eai + f3a 2 = 
/3a„_i - Ea„ + fia n +\ = 
(3a N -i - Ea N + 70^+1 = 
jaN - Ecin+i + PoaN+2 = 
/3oa„_i - Ea n + (3oa n +i = 



n<-l (CI) 
(C2) 
(C3) 

2 < n < N - 1 (C4) 
(C5) 
(C6) 

n > N + 1 (C7) 



a n being the coefficient of the wave function relative to 
the s atomic orbital centered on atom n. We look for a 
solution in which an electronic wave with crystal momen- 
tum fco > inside the left lead is reflected by the defect 
with a reflection amplitude probability r and transmitted 
into the right lead with a transmission amplitude proba- 
bility t. Then the amplitude of the electronic wavefunc- 
tion on each atom can be written: 




inko 
— ink 



if n < , 

if 1 < n < N , 

if n > N. 



(C8) 



which obviously satisfies Eqs IC1IC4I and IC7I with E = 
2/3oCOsfco = 2/3 cos k (Note that k may be imaginary 
when \P\ < |jSb|). Substituting [C8] for a n into Eqs E! 
IC3I IC5l and lC6l vields a linear system of 4 equations with 
the 4 unknown variables A, B, r and t. After solving this 
system the transmission factor is T^(E) = 2tt* given as 
a function of (3, (3q, 7 and N by Eq. [5] of the main text. 
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